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Abstract. The tropical arithmetic operations on R are defined by a © b = 
min{a, b} and a®fe = a + b. Let A be a tropical matrix and k a positive integer, 
the problem of Tropical Matrix Factorization (TMF) asks whether there exist 
tropical matrices B S R mxfc and C G R fexn satisfying B (g> C = A. We show 
that no algorithm for TMF is likely to work in polynomial time for every fixed 
k, thus resolving a problem proposed by Barvinok in 1993. TMF is also shown 
to be hard for matrices with bounded tropical rank. Proving that TMF can be 
solved by a polynomial-time algorithm if k < 3, we answer a question posed 
by Develin, Santos, and Sturmfels. Another question they have posed asks 
whether every tropical matrix of factor rank k has a rank-fc submatrix of size 
at most N(k) X N(k); we answer this question in the negative for every k > 4. 



1. Introduction 

The tropical semiring is the set R of real numbers equipped with the operations of 
tropical addition and tropical multiplication, which are defined by a(Bb = min{a, &}, 
a (g) b — a + b. The tropical semiring is essentially the same structure as the max- 
plus algebra, which is the set K with the operations of maximum and sum, and is 
being studied since the 1960's, when the applications in the optimization theory 
have been found [33] . The tropical arithmetic operations on R, which allow us to 
formulate a number of important non-linear problems in a linear-like way, arise in- 
deed in a variety of topics in pure and applied mathematics. The study of tropical 
mathematics has applications in operations research [T5] , discrete event systems [B] , 
automata theory [35] , optimal control [55[ 123 , algebraic geometry [TU [TB] , and oth- 
ers; we refer to [23] for a detailed survey of applications. A considerable number 
of important problems in tropical mathematics has a linear-algebraic nature. For 
instance, the concepts of eigenvalue and eigenvector, the theory of linear systems, 
and the algorithms for computing rank functions are useful for different applica- 
tions [3 HH [53] . Some applications also give rise to studying the multiplicative 
structure of tropical matrices [31], and in this context, the Burnside-type prob- 
lems are important [22j [32] . Another interesting problem is to study the subgroup 
structure of the semigroup of tropical matrices under multiplication [251 126] . 

In our paper, we consider the problem of matrix factorization, which is also 
related to the concept of factor rank of matrices over semirings 0. The study of 
factor rank dates back to the 1980's [9 , and has now numerous applications in 
different contexts of mathematics. Being considered on the semiring of nonnegative 
matrices, the factor rank is known as nonnegative rank and has applications in 
quantum mechanics, statistics, demography, and others [11] . The factor rank of 
matrices over the binary Boolean semiring is also called Boolean rank and has 
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applications in combinatorics and graph theory [7J [24] . Finally, for matrices over a 
field, the factor rank coincides with the classical rank function. 

In the context of matrices over the tropical semiring, the factor rank is also 
known as combinatorial rank [3] and Barvinok rank [15] , and the study of this 
notion has arisen from combinatorial optimization [2j. The factor rank appears 
in the formulation of a number of problems in optimization, for instance, in the 
Traveling Salesman problem with warehouses [4]. Also, the notion of factor rank is 
of interest in the study of tropical geometry [TJ1 [TS] , where the factor rank can be 
thought of as the minimum number of points whose tropical convex hull contains 
the columns of a matrix. Let us define the factor rank function for tropical matrices, 
assuming that multiplication of tropical matrices is understood as ordinary matrix 
multiplication with + and • replaced by the tropical operations © and (g>. 

Definition 1.1. The factor rank of a tropical matrix A £ R mx ™ is the smallest 
integer k for which there exist tropical matrices B £ R mxfe anc l C £ R fex ™ satisfying 
B <g> C = A. 

The most straightforward way of computing the factor rank is based on the 
quantifier elimination algorithm for the theory of reals with addition and order [17] . 
Indeed, Definition 11.11 allows us to define the set of all m-by-n matrices with factor 
rank k by a first-order order formula. We then employ the decision procedure based 
on the quantifier elimination algorithm provided in |17j to check whether a given 
m-by-n matrix indeed has factor rank k. 

However, the computational complexity of quantifier elimination makes the algo- 
rithm mentioned unacceptable for practical use. Another algorithm for computing 
the factor rank is given by Develin in the paper |14) , where he develops the theory 
of tropical secant varieties. He characterizes the factor rank from a point of view 
of tropical geometry, and the characterization obtained provides an algorithm for 
computing the factor rank. 

Unfortunately, neither the algorithm by Develin nor any other algorithm is likely 
to compute the factor rank of a tropical matrix in polynomial time. Indeed, the 
problem of computing the factor rank is NP-hard even for tropical 01-matrices, 
see |15j . In other words, the general problem of Tropical Matrix Factorization 
( TMF), which asks whether a given matrix A and a given integer k are such that 
A = B ® C for some B £ R mxfe and C £ R kxn , turns out to be NP-hard. 

Besides the general problem of computing the factor rank, the problem of detect- 
ing matrices with fixed factor rank deserved attention. Certain hard problems of 
combinatorial optimization admit fast solutions if the input matrices are required 
to have factor rank bounded by a fixed number [2] [3]. The Traveling Salesman 
problem (TSP) also admits a fast solution if we require the distance matrix to 
have a fixed factor rank [?J [5] ; this special case of TSP is also known as TSP with 
warehouses [4] . Matrices with bounded factor rank arise naturally in the problems 
mentioned and in a number of other problems in combinatorial optimization [2] 
and tropical geometry [15]. These considerations led to the following interesting 
question on tropical matrix factorizations. 

Question 1.2. [H|3] Does there exist an algorithm that solves TMF for every fixed 
k in polynomial time? 

In [2 , Barvinok expected that Question [O] can be answered in the positive, this 
question has also been mentioned in [3]. Further investigations on the problem of 
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determining matrices with fixed factor rank have been carried out in \10\ I13| If 4[ 
[HI [20]. However, the problem remained open even in the case k = 3 and has been 
formulated again in |15) . 

Question 1.3. [HI Section 8, Question 2] Is there a polynomial-time algorithm 
that decides whether a given tropical matrix has factor rank 3? 

In [15], a number of related problems has also been posed. The following question 
asks whether the set of matrices with fixed factor rank can be defined by ranks of 
the submatrices of bounded size. 

Question 1.4. [T5J Section 8, Question 3a] Is there a number N(k) such that if 
all minors of a tropical matrix A of size at most N(k) have factor rank at most k 
then A itself has factor rank at most kl 

Assume that Question 11.41 is answered in the positive for some k. To solve 
Problem 12. 1[ it is then sufficient to check whether every minor of size at most 
N(k) x N(k) has factor rank at most k. The number of minors of fixed size is 
easily seen to be polynomial in the size of a matrix, and each of these minors can 
be treated in fixed time thanks to, say, quantifier elimination algorithm. Thus we 
see that, indeed, as noted in [15] . a positive answer to Question 11.41 would allow to 
construct a fast algorithm for recognizing tropical matrices with factor rank k. 

Besides the factor rank, there are other important rank concepts for tropical 
matrices, see [TJ [15]. One of them is presented by the function of tropical rank, 
which arises naturally in tropical linear algebra and tropical geometry |I5| . First, 
we define the tropical permanent (which is also called the tropical determinant) of 
a tropical matrix D £ M. nxn as 

(I.I) perm(D) = min {£>i,<t(i) + • • • + D n . a ( n) ) , 

where S n stands for the symmetric group on {1, . . . ,n}. D is said to be tropically 
singular if the minimum in is attained at least twice. The tropical rank of 

a tropical matrix A is the largest integer r such that A has a non-singular r-by-r 
submatrix. Now we can formulate another problem related to the complexity of 
tropical matrix factorization. 

Question 1.5. [HI Section 8, Question 3b] Is there a polynomial-time algorithm 
for the factor rank of matrices with bounded tropical rank? 

The progress in solving the problems we have mentioned has mostly been based 
on studying matrices with factor rank at most 2, and the set of such matrices is now 
indeed well studied. First, the TMF problem with k < 2 can be solved by a linear- 
time algorithm, see (TO] [15] . Further, it has been proven in [15] that the factor rank 
of a matrix is at most 2 if and only if all its 3-by-3 minors have factor ranks at most 
2. The set of d-by-n matrices with factor rank 2 has been studied as a simplicial 
complex in [13]. This set has also been studied from the topological point of view 
in [2D], and the space of d-by-n matrices of factor rank two modulo translation 
and rescaling has been shown to form a manifold. For k > 3, the question of fast 
algorithm for recognizing tropical matrices with factor rank k remained open. 

In our paper, we answer Question 11.21 showing that no algorithm is likely to 
solve TMF in polynomial time for every fixed k. More precisely, we show that the 
problem of recognizing tropical matrices with factor rank 8 is NP-hard. However, it 
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turns out that the answer for Question ll.3l is positive. We also answer Question QT4] 
in the negative for every k > 4 and obtain a negative answer for Question ll.51 

The following notation is used throughout our paper. By Ay we will denote the 
(i, j)th entry of a matrix A, by A[r±, . . . ,r p \c±, . . . ,c q ] the submatrix of A formed 
by the rows with labels r\ , . . . , r p and columns with c\ , . . . , c q . We will say that 
tropical matrices B,C G s. mxn coincide modulo scaling if there exist real numbers 
a%, . . . , a m , /3%,.. . ,0 n such that Hy = Cy + cti + j3j for all i and j. Straightfor- 
wardly, B and C have the same tropical and the same factor ranks in that case. 

2. Preliminaries 

In this section, we give a more specific formulation for the TMF problem and 
prove the preliminary results on computational complexity of this problem. We 
describe the construction of the extended tropical semiring KU {oo}, which will be 
helpful for further considerations, and we prove a number of auxiliary results on 
the factor rank of matrices over R U {oo}. 

2.1. Tropical matrix fc-factorization. We will consider a family of decision prob- 
lems which specify a value of inner dimension k in the general formulation of the 
TMF problem. We require an input matrix to consist of integers in these prob- 
lems, since this assumption will be convenient in our further considerations. For a 
positive integer k, the problem we deal with is as follows. 

Problem 2.1. TROPICAL MATRIX fc-FACTORIZATION (fc-TMF). 
Given a tropical matrix A G Z mx ™. 

Question: Do there exist tropical matrices B G R mxfe anc j Q G R fcXTl satisfying 
B ® C = A? 

In the input of Problem 12.11 as well as throughout our paper, we assume that 
an integer matrix is presented in the usual way as a rectangular array of numbers 
written in the decimal system. Namely, matrices will be considered to be words over 
the alphabet that consists of the characters "0", "1", . . ., "9", "+", a space, 
and an cnd-of-line symbol. In particular, the size of the input in Problem 12 . 1 1 equals 
Y^, [log io rnax{|Ay |, 1}] + 3mn, where [x] denotes the integer part of x. 

Now we will prove that the TMF problem belongs to the class NP if an input 
matrix is required to consist of integers. We need the following technical lemma. 

Lemma 2.2. Let a matrix A G Z mxn have factor rank at most k. Let g and I 
denote, respectively, the greatest and the least elements of A, set also h = \g\ + |/|. 
Then there exist matrices B G Z rnxk and C G Z kxn such that B ® C — A and 
\Bi T \ < h, \C T j \ < h for any triple of indexes (i,j,r). 

Proof. By Definition 1 1.11 there are matrices B' G R mxk and C G R kxn satisfying 
B' ® C = A. Consider matrices B" and C" defined by B'( T = B[ T - {B' lT }, 
C"j = C' T j + {—C' T j}, where {x} = x — [x] is the fractional part of x. The matrices 
B" andC" are then integer, and we have \(B? T + C^) - (B' lT + C' T:j )\ < 1. Therefore, 
B' iT + C' TJ = Aij implies B" r + = A tJ , and from B' iT + G' Tj > Mj it follows that 
B" T + C'^ > A^. Thus we have that B' ® C = A implies B" ® C" = A. 

Further, we subtract b T , the least element of the rth column of B", from every 
entry of the rth column of B" and add b T to every entry of the rth row of C". The 
matrices B and C obtained satisfy B ® C = A, and zero appears as the minimal 



THE COMPLEXITY OF TROPICAL MATRIX FACTORIZATION 



5 



element of every column of B. The definition of matrix multiplication then shows 
that every entry of C is greater than or equal to I. Finally, those entries of B and 
C that are greater than \g\ + \l\ can be then replaced by \g\ + \l\ without changing 
the product B ®C. □ 

Now we can prove our first results concerning the computational complexity of 
tropical matrix factorization. 

Theorem 2.3. Given a tropical matrix A £ 2 mxn and a positive integer k. The 
problem of deciding whether A has factor rank at most k belongs to the class NP. 

Proof. The factor rank of A is at most min{r7i, n} (see [HI Proposition 2.1]), so the 
problem can be solved immediately if k > min{m, n}. For k < min{m, n}, we apply 
Lemma \2. 21 If A has factor rank at most k, then there exist matrices B £ z mxfc 
and C £ Z' cxn that satisfy B ® C = A and require at most (m + n)k (log 10 h + 3) 
symbols to be written (h denotes here the sum of the absolute values of the greatest 
and the least elements of A). The size of the input data of the problem considered 
is at least log 10 -| + mn, so the result follows. □ 

Theorem 12.31 shows that TMF, the general problem of tropical matrix factoriza- 
tion, belongs to NP if the input matrices are assumed to consist of integers. The 
problem fc-TMF, introduced in this section, is thus in NP as well. 

Theorem 2.4. The k-TMF problem belongs to NP. 

Proof. Follows from Theorem 12.31 □ 

2.2. Extended tropical semiring. It will be useful for our further considerations 
to extend the tropical semiring by an infinite positive element, which we denote by 
oo. We also write R for R U {oo} and assume a © oo = a, a ® oo = oo for a £ R. 

Defining the factor rank, one can now think of a tropical matrix A £ R mxn as a 
matrix over R, and allow matrices B and C from Definition 11.11 to contain infinite 
elements. However, it turns out that the rank function defined in this way is the 
same as that defined with respect to Definition ll.il Indeed, if matrices A £ R mxn , 
B £ R , C £ R satisfy A = B <S> C , then we can replace infinite entries of B 
and C with a sufficiently large real without changing the product of the matrices. 

We can therefore extend Definition 11.11 to the case of matrices over R and define 
the factor rank of A £ R to be the smallest integer k for which there exist 
B £ R mX and C £ R satisfying A — B ® C. The following technical lemma 
will play an important role in the considerations of our paper. 

Lemma 2.5. Let zero appear as a minimal element of every row and every column 
of a matrix A £ R . Consider the matrix A' obtained from A by replacing every 
infinite entry with the number 2g + 1, where g stands for the maximal finite entry 
of A. Then the factor ranks of A and A' are the same. 

Proof. Consider matrices B' £ R mxfc and C £ R fcx ™ satisfying B' <g> C = A'. For 
t £ {l,...,fc}, we subtract bt, the least element of the tth column of B, from every 
entry of that column and add bt to every entry of the tth row of C . The matrices 
B" and C" obtained satisfy B" ® C" = A' , and zero appears as a minimal element 
of every row of B" and of every column of C" . Further, we replace by oo every 
entry of B" and C" that is greater than g, and we denote the matrices obtained by 
B and C. It is then easy to check that B (g> C = A. 
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So we have proven that the factor rank of A is at most that of A' . It is thus 
sufficient to note that A' = E ® A, where E stands for an m-by-m tropical matrix 
satisfying Eij = if i = j and Efj = 2g + 1 otherwise. □ 

The following easy observation will also be useful. 

Proposition 2.6. The factor rank of a tropical matrix A E W mxn is one less than 
that of the matrix A' E R(m+i)x(n+l) defined by A^ = A ij} A' l n+l = A' m+1J = oo 
for i<m, j <n, and A' m+l n+1 = 0. 

The rest of our paper is organized as follows. In Section 3, we consider the fc- 
TMF problem for large enough values of inner dimension fc. We show that fc-TMF 
is NP-complete for k > 8, thus giving a negative answer for Question 11.21 As a 
corollary of this result, we show that the problem of deciding whether the factor 
rank of a matrix A exceeds 8 is NP-hard even if A is assumed to have tropical rank 
at most 8; in particular, this gives a negative answer for Question ll.5l Note that this 
situation is in contrast with that of the Gondran-Minoux rank, which can be com- 
puted in polynomial time if a given matrix has bounded tropical rank [18] . In Sec- 
tion 4, we consider the problem of tropical matrix factorization for small values of k. 
We construct an algorithm that requires at most O ((m + n) 3 mnlog(mn) + m 3 n 3 ) 
arithmetic operations and decides whether a given matrix A 6 R rax " has factor- 
rank at most 3. This answers Question 11.31 in the affirmative and shows that the 
3-TMF problem belongs to the class P. In Section 5, we consider the question of 
whether matrices with fixed factor rank are defined by ranks of their bounded-size 
minors, and we give a negative answer for Question II .41 in the case when k > 4. 

3. 8-TMF IS NP-COMPLETE 

In this section, we consider the problem of tropical matrix factorization for suf- 
ficiently large but fixed values of fc, and show that this problem remains hard if 
the dimension k is assumed to be bounded. We show that no algorithm is likely to 
solve fc-TMF efficiently if k > 8. More precisely, we construct a polynomial-time 
reduction from an NP-complete problem to fc-TMF, thus proving that fc-TMF is 
NP-hard. In other words, we show that efficient solution of fc-TMF would provide 
efficient solutions to many problems that are of interest for a long time and widely 
believed to be hard [21] . Therefore, the results of this section give a negative answer 
for Question II. 21 As a corollary, we obtain a negative answer for Question II. 51 

The NP-complete problem which we reduce to fc-TMF is known as SET SPLIT- 
TING. This problem has been considered in [21] and is formulated as follows. 

Problem 3.1. SET SPLITTING. 

Given a set S and a family C of subsets of S. 

Question: Are there sets <&i and <& 2 such that $i R $2 = 0, $1 U $2 = S, and 
no subset from C is entirely contained in either $1 or $2? 

3.1. Reformulation of SET SPLITTING and some specific conventions. 

Our reduction to fc-TMF, however, does not go from SET SPLITTING directly. 
Let us introduce the following reformulation of SET SPLITTING, which is useful 
for our considerations. 

Problem 3.2. SET SPLITTING REFORMULATED. 

Given a nondecreasing integer sequence gq = 0,eri, . . . ,er m and a sequence of 
pairwise disjoint sets . . . , lZ n such that IZi U . . . U lZ n — {!,..., cr m }. 
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Question: Are there sets $1 and $2 such that $1 n$ 2 = 0, $1 U<1>2 = {1, • • • , n], 
and the sets 



are both nonempty for every /i G {1, . . . , m}? 

Problem EH as well as SET SPLITTING, turns out to be NP-complete. 
Lemma 3.3. Problem \3.S\ is NP-complete. 

Proof. Let S = {si, . . . , s n } be a set and Ci, . . . ,C m a sequence of its subsets, by 
iff we denote the set of all pairs (si,Cj) satisfying s,; G Cj. We set ao = and 
°V = \Ci \ + ••• + \Cfj,\ for every \x G {l,...,m}. Then it is easy to construct 
a bijection v : iff — > {l,...,er m } mapping every (si,Cj) G iff to some integer x 
satisfying <7j_i < x < Uj. Finally, for i G {1, . . . , n}, we define IZi to consist of all 
numbers whose preimage under v has as the first entry. These settings can be 
straightforwardly checked to give a polynomial reduction from SET SPLITTING 
(which is known to be NP-complete, see [21]) to Problem 13.21 On the other hand, 
the fact that Problem 13.21 is in NP is clear. □ 

Our way of proving that 8-TMF is NP-hard is based on the polynomial-time 
reduction from Problem [32] In order to describe this reduction explicitly, we need 
to specify some conventions that will be used throughout this section. We will 
denote the input data of Problem l3.2l bv E = (ao, . . . , a m ), 1Z = (Hi, ■ ■ . , H n ). The 
conventions to be used are given by the following remarks. 

Remark 3.4. Note that if |<7j — o~j\ < 2 for some distinct i and j, then the answer 
for Problem 13.21 is trivially 'no'. Consequently, Problem GO] remains NP-complete 
under an additional assumption that \o~i — o~j\ < 2 implies i = j. In what follows, 
we suppose that S satisfies this assumption. 

Remark 3.5. The formulation of Problem 13.21 requires that for every u € 
{1, . . . , er m }, there exists a unique v — v(u) such that u G 1Z V . 

Remark 3.6. For any u G {1, . . . , <r} and v G {1, . . . , a — m}, we set 



where A and (i stand for the maximal subscripts satisfying u > <t\ and v > — /i, 
respectively. We will also write cr for a m , and G for the number 21a. 

3.2. Encoding the input of Problem 13.21 by a matrix. Let us now begin 
to describe the reduction from Problem 13.21 to 8-TMF. Namely, this subsection is 
devoted to the construction of a tropical matrix A — A(S, 1Z) which we assign 
to the input data (E,1Z) of Problem 13.21 By our setting, the matrix A will have 
cr — m + n + 9 rows and a + 10 columns. The following convention will be useful to 
present this matrix in a more concise way: The rows of A will be labeled with the 
elements of the set 




(3.1) 



= 20u - 10A - 9, p. 



20v+ 10/x + l, 



I = {!♦, . . . , 94, 10, ■ • ■ , n<>, IV, . . . , (a - m)V}, 



the columns with the elements of 



J = {l^...,10Mtf,...,crfl}. 

Now we can define the entries of ^4(S, TV). 
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Definition 3.7. Let (£,72) be the input of Problem \3.2\ u be an arbitrary number 
from {1, . . . , er}, v from {1, . . . , a — to}, r] from {1, . . . , n}. By p u , j v , and G we 

denote the numbers from Remark \3.b\ The entries of the matrix A(H,TV) G R 
are defined as follows. 

(1) If u and v are such that u = v + p = for some p G {1, . . . , m}, then 
A v y lU i — min{7„, p v , 10(u + v) — 11}. Otherwise, if no p satisfying that 
condition exists, then A v ^>^ — min{7„, p v , 10(u + v) — 10}. 

(2) A v $. u $ = min {7^, G — 2G \rj — v(u)\} (recall that v{u) denotes a number 
such that the condition u G 72, y ( M ) is satisfied). 

(3) Ayvj^ — 0, A v v^ — p v , A v y^ = Aav^ = 00; A v v tT ^ — 2 for r G 
{1,...,6}. 

(4) = = A v $,th = 2, Aj^.s^ = -2nG, A^^ = 2nG, A^^ = 
2(n + 1 - r])G; A v<> ^ = for r' G {2, 3, 4, 5}. 

(5) ^4i<|»,u)j=^2*,«tl = ^3*,utt=2, At*,uj = 0, ^6*,utl = 10w - 3, A74 itt j = 7„, 
^8*,uij = -2Gi/(«) + G, A 94lkM = -2G(ri + 1 - + G. 

(6) If u = cr^-i + 1 for some p G {1, . . . , to}, i/ien As^^j = 10u + 1; otherwise, 
^5#,«)t = 10it - 3. 

(7) TTie matrix A[14t, . . . , 94k|lt], . . . , 101]] equals 



(3.2) 



/ 


2 





2 
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00 


00 


\ 
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00 


00 
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00 


00 
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00 
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00 


00 


00 






00 
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00 


00 


00 


00 


00 


00 





00 




V 


00 


00 


00 


00 


00 


00 


00 


G' 


00 





/ 



where G' stands for —2(n + 1)G. 

It is easy enough to show that the function assigning the matrix A(£,72) to a 
tuple (£,72) has the polynomial complexity. 

Lemma 3.8. Given a tuple (£,72). The number of operations needed to construct 
the matrix A(£,72) is polynomial in the size of the data. 

Proof. The numbers of the rows and columns of A(£,72) are given by linear func- 
tions of variables appearing in the data. On the other hand, every entry of A(£, 72) 
can clearly be computed in bounded number of operations. □ 

We will further consider tropical factorizations of the matrix A = A(£,72). In 
order to adapt the notation of the row and column indexes of A to the product of 
matrices, we need the following convention. 

Remark 3.9. Throughout this section, for considering matrices B and C satisfying 
A(£, 72) = B ® G, we will use the following notation. The elements of X will 
label the rows of B, the elements of J the columns of G, the numbers 1, . . . , k the 
columns of B as well as the rows of G. In our settings, A — B <g> G means that 
Aij = mm^. =1 {Bi T + C T j} for all i G I and j G J . 

The following lemma shows that the factor rank of A(£, 1Z) does not exceed 8 if 
the question of Problem 13.21 is answered in the affirmative on the query (£,72.). 
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Lemma 3.10. Let an input (S,TZ) be such that the answer in Problem \3.2\ is 'yes'. 

|X|x8 8x|J7" 

Then there are matrices B £ R and C € K such that A(T, 1 TV) — B ® C . 

Proof. The assumption of the lemma implies that there exist $1 and $2 such that 
$1 n $ 2 = 0, $1 U $ 2 = {1, • ■ • , n}, and the sets 



U n i f|W-i + 1. • • • . ^1 and |J UA f|{^_ x + 1, 



J 



Vje*! / \ie c i > 2 / 

are both nonempty for every p, 6 {1, . . . , m}. 

Then for every p, the set -{V(<7 M _i + 1), . . . , ^(c^)} has nonempty intersections 
with both $1 and $2 (recall that, with respect to Remark 1331 v{u) denotes a unique 
number for which the condition u 6 TZ V ( U ) is satisfied). By "H^ we will denote the 
greatest number of the set {ct^-i + 1, . . . , for which the numbers v^M.^) and 
vicrfj) are neither both in $1 nor both in $2- 

Now we can present the matrices B and G. In the rest of our proof, we will 
denote an arbitrary number from {1, . . . , <r} by u, from {1, . . . , a — m} by v, from 
{1, . . . , n} by 77; also, \ wml stand for either 1 or 2. As in Definition 13.71 p Ul ~{ v , 
and G stand for the numbers defined in Remark 13.61 In items B4, B5 and C3, we 
assume that the index p runs over the set {1, . . . , m}. 



Bl: The matrix S[l*,...,9*|l, 



is equal to 



(3.3) 



V 



2 


2 





2 


2 


2 


CO 


CO 





2 


2 


2 


2 


2 


CO 


CO 


2 





2 


2 


2 


2 


CO 


CO 


2 


2 


00 





2 


2 


CO 


CO 


00 


00 


00 


CO 





2 


CO 


CO 








00 


CO 








CO 


CO 


00 


00 


00 


CO 


CO 





CO 


CO 


00 


00 


00 


CO 


CO 


CO 





CO 


00 


00 


00 


CO 


CO 


CO 


CO 






if 77 € <f> 



X' 
= 0°, 



and B v<> , x 
= 2, 



= co if 77 
B,,0,7 = 2tjG, and ^0,8 



<& x ; take also 



B2: We take B v<>tX 
2(n + l-r?)G. 

B3: We set B v y >3 = B v y i7 = B v y jS = 00, B v v A = B v ^r = p v , B v <y j6 = 0. 
B4: If g G {07,-1 — /i + 2, . . . , <t m — /i} is less than — fi + 1, then we set 

B g v tX = 10(7 — 7; if 5 is greater than — /i + 1, then B g <y^ x = IO5 — 8. 
B5: Finally, if g = — p, + 1, then we take B g ^. x = 10g — 8 if v{T-L^) £ & x , 

and S 99 , x =' 10a - 7 if i $ x . 

We proceed with the description of the matrix G. 



(3.4) 



matrix C[l, . 


..,8 




. . , 10lq] is 


equal to 










/ 


2 








2 


2 


2 


2 


00 


00 


00 


\ 






2 


2 





2 





2 


2 


00 


00 


00 








2 





2 








2 


2 


00 


00 


00 











2 


2 


2 


2 


2 


2 





00 


00 








2 


2 


2 


2 


2 





2 





00 


00 








2 


2 


2 


2 


2 


2 





00 


00 


00 








00 


00 


CO 


CO 


CO 


CO 


CO 


00 





00 






V 


00 


00 


CO 


CO 


CO 


CO 


CO 


G" 


00 





/ 
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with G' standing for — 2(n + 1)G. 
C2: We set C^j = C6.«j = 7«, Giutf = 0, C5.„j = As*,^, CV.^fl = -2v{u)G- 
G, C 8 , ua = -2(n + 1 - v(u))G + G. 



C3: Assume h G {o>-i + 1, . . . , cr M }. We set C 



x-H 



if v{h) G 5? x , set also 



C 



10ft. - 2 if v(h) £ <I> X and < h < a^, finally, set C x ^ = 10ft - 3 



if v(h) <E> X and either ft < or ft = 



It is straightforward to check that the matrices B g 1 and C G 
are well defined. In order to complete our proof, we need to verify that indeed 
B®C = A(E,n), that is, 



(3.5) 



At 



min{i? 

T — 1 



+ C Tj } 



for alii G I and j G J . Consider the six possible cases. 

1. For i G . . . , 94*} and j G {It], . . . , 10t]}, it suffices to check that the product 
of the matrices (|3.3j) and (|3.4j) equals the matrix f|3 . 2j) . 

2. Let i = 7/0 and j G {It], . . . , 10t]}. We need to verify that 
the row (A v q^, . . . , A v $ t io§) is a product of the row of B labeled by f?<> 
and the tropical matrix (|3.4p . Namely, we need to check that the row 
(B^i, Btjo.2, 0, oo, oo, 2, 2f?G, 2(n + 1 — ?y)G) multiplied by the tropical ma- 
trix fliS} equals (2, 0, 0, 0, 0, 2, 2, -2r]G, 2rjG, 2(n + 1 - rj)G). We note that also 
{B v Q t i, B n ^ i2 } = {0, oo} by the definition of B, checking is then straightforward. 

3. For i = •y's? and j G {lt],...,10t]}, checking (|3.5[) means checking that the row 
(A v y t i^ . . . , At^ioii) i s a product of the row of B labeled by and the tropical 
matrix (|3.4[) . We note that miri{ i, -B„<?. 2 } > 2 by the definition of B. Then 
the row (B v <? t i, B v y t 2, 00, p u , p„, 0, 00, 00) multiplied by the tropical matrix (|3.4j) 
can be straightforwardly checked to equal the row (2, 2, 2, 2, 2, 2, 0, /0„, 00, 00). 

4. Further, assume i G {l<|k, . . . , 9d(k} and j = lift. Observe that we have 
min{Gi, U |j, G2,n|j} = 10ti — 3 for it = 07,-1 + 1, and A^^ = lOu — 3 for 
u 7^ cr M _i + 1. Also, by our settings, Gi )U j, C2,uj) G {10u — 2, 10w — 3,7„} and 

5*,u|j G {10u — 3, lOu + 1} for any u. 



A 



So we can check that 





2 


2 





2 


2 


2 


oc 


oc 


\ 







2 


2 


2 


2 


2 


oc 


oc 






2 





2 


2 


2 


2 


oc 


00 






2 


2 


00 





2 


2 


00 


00 






00 


00 


00 


00 





2 


00 


00 












00 


00 








OO 


00 






00 


00 


00 


00 


00 





00 


00 






00 


00 


00 


00 


00 


00 





00 




V 


00 


00 


00 


00 


00 


00 


00 





/ 



® 



V 



Gi,„j 

In 



A5*,«)l 

7« 
Gi 

G 2 



V 



2 
2 
2 


10m -3 

lu 

G'i 
G 2 



\ 



where Gi and G 2 stand for —2Gv{u) + G and —2G(n + 1 — f (u)) + G, respectively. 
This implies that (|3.5[) holds. 

5. Now consider i = t/<0> and j = ujj. In this case, we have min{i?,^.i + 
Ci.uftj^O.a + C 2m} > °) and min{S r) ^,i + Ci iUt t,B,o )2 + C^utf} = 7„ if ?? = 
We also have min{_B r; ^> ! 7 + Cj^, B v ^ t s + Cs.uj)} = G — 2G |?/ — and 
Br]fy,T' + Ct'.uJ > 7u f°r t' G {3,4,5,6}. This implies 



min{B r) ^. T + G r . utt } = min{7„, G - 2G \q - u(u)\} . 

T = l 
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thus proving p.5[) . 

6. Finally, we assume i = v^?, j — ujj. In this case, the minimum in Q3.5P 
computed for r running over {3, . . . , 8} equals xma{p v , 7„}. It is therefore sufficient 
to prove that 

In] 

for all u and v. Note that by the definitions of B and C, it holds that B v <y a +C XtU $ > 
10m + 10f — II for x G {1: 2}. By A and \i we will denote the maximal subscripts 
satisfying u > <J\ and v > — /i, respectively. Let us consider possible cases. 

6.1. If lOu + lOv — 11 > min{p„,7 u }, then the minimum in (|3.6p equals 
min{p„,7 u }. Item 1 of Definition 13.71 shows that also j4„<?,utt = min{p„,7 u }, prov- 
ing 4SSD. 

6.2. If the assumption of item 6.1 fails to hold, then lOu+lOv — 11 < 20u— 10A — 9 
and 10it + 10u-ll < 20v + 10// + 1, so that 10t; + 10A < 10m + 2 < 10u+ 10/i+14, 
and we obtain v + /i + 1 > u > v + A. In particular, we see that the assumption of 
item 6.1 holds if A > /i + 1, so (|3 . 6[) holds in this case as well. 

6.3. Now suppose A < /i. The definition of A and p implies u < ax+i — 1 and 
v > ex n — /i. By Remark 13.41 — ux+i > 2(/i — A — 1). Thus we have 

U < (Tx+l — 1 < <T\+1 + fl-\-2<(T^-fI+X<V + X, 

and then lOv + 10A > lOu + 2, so lOw + lOv — 11 > j u . We see that the assumption 
of item 6.1 holds again, proving (|3.6p . 

6.4. Further, suppose A = /i + 1. By items 6.1 and 6.2, it can be assumed without 
a loss of generality that u = v + A. In this case, the definition of A and p shows that 
u> a\ and v < a\ — A, so we have u = v + A = a\. Note that v{%\) G $ x if and 
only if v{(J\) ^ $ x by the definition of the index H\. From items B4, B5 and C3 
it therefore follows that min{B v ^,i + Ci.uj, B v <y^ + C2,«|)} = lOti + lOt; — 11, from 
item 1 of Definition 13.71 that A v <y tU $ — min{10(u + w) — ll,j u ,p v }, so (|3.6p holds. 

6.5. Finally, assume X — p. We use items 6.1 and 6.2 and then assume without a 
loss of generality that either u = v + fJ,oiu = v + fi + l. In this case, the definition 
of A and p shows that u G {u^ + 1, . . . , c^+i — 1}. Now from items B4, B5 and C3 
it follows that minl-B^i + Ci, M j, B v <$^ + C2,utt} — 10u + 10u — 10, from item 1 of 
Definition 13.71 that A v <y^ — min{10(M + v) — 10, -f u ,p v }, so (|3.6[) holds again. 

Items 6.1-6.5 consider all the possibilities, so the consideration of case 6 is com- 
plete. The cases 1-6 prove (|3.5p for every (i,j) £ I x J, so we can conclude that 
indeed B ® C = A{H,1l). □ 

3.3. The main results. In this subsection, we finalize the proof of the fact that 
fc-TMF is NP-hard for k > 8. What we need is the following lemma. 

Lemma 3.11. Let an input (T,,1Z) be such that the answer in Problem \3.S\ is 'no', 
assume B eR mx& , C &R SxlJl . Then B ® C ^ A(T,,TZ). 

Proof. The proof is by reductio ad absurdum. We assume that B ® C = A(E,1Z), 
in other words, that the equality 

8 

(3.7) Aij = min{B 4T + C rj } 

T—l 

holds for all i G I and j G J. 
Step 1. 
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1.1. Let t' £ {1, ...,8} satisfy C T >^ ^ oo for some tp £ {It], . . . , 7t|}. From 
Definition 13.71 it follows that As4t,ip = -<4-9*,V> = oo, so the equality (|3.7|) implies 
B&4,, T ' + C t >tp = Bg^, iT r + C T 'j(, = oo, and thus we have i?s*,r' = -39*, t' = oo. 

1.2. By t — {t%, . . . ,T g } we denote the set of all r £ {1, ... ,8} satisfying both 
min{C T) 9|], C Ti io^} ^ oo and C T .\^ — . . . — C T ,7t = oo. From item 1.1 it follows that 
either -Bs*,t" = Bg^ tT n — oo or C T " ,% — C T ",ioti ~ 00 f° r ever y r " ^ t. Thus the 
equality (J377J) implies that 

A[84, 94|9^, 101]] = B[8*,9*|ri, . . . ,T g ] ® C[t x , . . . ,T g \% 101;], 

and, since from Definition 13.71 it follows that A[84, 94|9t], 10tj] = < J > ), we see 
that t has cardinality at least 2. 

1.3. Using item 1.2, we assume without a loss of generality that 7,8 £ f. Thus, 

(3.8) min{C 7 ,9i,, C 7i iot]} ^ oo, min{C 8 , 9 t,, C 8 ,ioti} °o, 

(3.9) C 7i , = C 8 i, = oo for ^ G . . . , 7^}. 
Step 2. From the equality (|3.7[) it follows that the condition 

8 

(3.10) mm{A s ,^, A s > w $ = min {B s /, r + mm{CV i9t] , CVioJ} 

r— 1 

holds, in particular, for every s' £ {14, • • • , 74, l'v', . . . , (cr — m)9}. 

By Definition 13.71 we have A s i gh — ^4 s ',iot] = °°> so that the left-hand side 
of (|3.10[) is indeed equal to oo. Now from the conditions (|3.8[) it follows that 

(3.11) B s , 7 = B y8 = oo for s' £ {1*, . . . , 7*, 1^, . . . , (<r - m)^}. 

Step 5. Note that the product B <E> C will remain unchanged if we add some 
number rgRto every entry of the rth column of B and —r to every entry of the 
Tth row of C. So we can assume without a loss of generality that for t £ {1, . . . , 6}, 
the tth row of C either consists of oo's or has as its minimal element. 

3.1. If the tth row of C consists of oo's, then, analogously, we can add a suffi- 
ciently large number to every entry of the tth column of B without changing the 
product B (g) C. Thus it can be assumed without a loss of generality that in this 
case, the tth column of B consists of positive elements. 

3.2. Now assume that the tth row of C has as its minimal element (namely, 
C tg ' = for some g' £ J) and that B s ' t < for some s' £ {14, . . . , 74, 19, . . . , (cr — 
771)9}. Then from the equality (|3.7|) it follows that A s i g i < B s * t + C tg ' < 0, a 
contradiction with Definition 13.71 

As a conclusion of Step 3, we therefore obtain that B s i t > 0, Ctj > for all 
s' £ {14, . . . , 74, 19, . . . , (cr - m)9}, t £ {1, . . . , 6}, and j £ J. 

Step 4 • Fix arbitrary u £ {1, . . . , cr} and v £ { 1, . . . , cr — m}. 

The condition (|3.1ip implies that mmf =1 {_B,i< t + Ctj} = min® =1 {B s r T + C T j} for 
s' £ {14, ■ • ■ , 74, 19, . . . , (cr — m)9}. In particular, this means that the product of 
matrices B = £?[14, • • • , 74, v9|l, . . . , 6] and C = C[l, . . . , 6|lt], . . . , 8t|, uft] is equal 
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to the matrix A[l4k, ■ ■ ■ , 74, v^\l\\ 



By Definition 13.71 we have that 



(3.12) 



B®C = 



t 


* 













* 


* 


* 


* 


\ 




* 








* 


* 


* 


* 


* 


* 






* 


* 





* 





* 


* 


* 


* 









* 




* 


* 


* 


* 












* 


* 




* 


* 





* 





* 






* 






















* 






* 


* 


* 




* 


* 





* 


* 




\ 


* 


* 


* 


* 




* 





* 


* 


/ 



where the asterisks stand for certain (possibly different) positive elements from K. 
By the conclusion of Step 3, the matrices B and C consist of nonnegative elements. 
Then, it turns out that the equality f|3 . 1 2[) determines the 'zero patterns' of B and 
C uniquely up to a permutation of row and column indexes. Namely, one can check 
that 



(3.13) 



fl[l*,...,7*,»<?|w(l),...,w(6)] = 



* 
* 
* 


o J 



(3.14) C[w(l),...,w(6)|lll,...,8l],«|t] = 



* 




* 




where the asterisks stand again for certain positive elements from R, and uj for some 
permutation on {1, . . . , 6}. Note that performing the same permutation of column 
indexes of B and row indexes of C will not change the product B ® C. Thus we 
can assume without a loss of generality that w is the identity permutation. 

Step 5. Fix arbitrary rj e {1, . . . , n} and h 6 {2, 3, 4, 5}. 

5.1. From Definition 13.71 it follows that = 0. 

5.2. The condition (|3.9p shows that CV^ = Cg^ — oo. 

Further, assume that -Bj^.t' + Cf./ifo < for some t' s {4, 5, 6}; then from (13. 14)) 
it follows that B v ^ t > < 0. In this case, we use (|3.14[) again and obtain that 
B V<>A + C^ X \ < 0, or B v<> ^ + C^et, < 0, or else B v<>fi + C 6Jii < 0. On the 
other hand, Definition 13.71 implies A^^^ — A^^ = ^0,7^ = 2, so we get a 
contradiction with the equality (|3.7|) . Thus we have m.m t „ =4 {B v $ :t " + CV'jit]} > 0. 

5.3. Items 5.1 and 5.2 and the equality (|3.7j) show that miri( =1 {B r) ^ it +Cf j^} = 0. 
From the equality (|3.14|) it then follows that 



(3.15) B v<><3 = 

and mm{B v<>tl , £^0,2} = 0. 



for rj e {1, 



'}• 



14 



YAROSLAV SHITOV 



5.4. By "fi (respectively, ^2) we denote the set of all r\ satisfying i?^o,i = 
(respectively, B^.2 = 0); from item 5.3 it follows that $iU$2 = {lv, n}. The 
initial assumption of the lemma states that, given the input (£,7£), the answer in 
Problem l3.2l is 'no'. Therefore, for some /t € {0, . . . , m — 1} and x 6 {1,2}, we have 

(3.16) {a M + l,...,<r M+ i} C (J Hi- 

5.5. The condition (|3.16j) means that for every u' G {er M + 1, . . . , cr^+i}, there 
exists r)(u') G {1, . . . , n} for which u' G Tl v (u') an d v( u ') 6 ^x- (He-re and in the rest 
of the proof, /i and \ denote the indexes defined in item 5.4.) Item 2 of Definition l3.7l 
then shows that A^^i^^^ = , the definition of ^ x that B v ^,^ jX = 0. The 
number j u / is to be understood here with respect to (|3.1[) . namely, for our case, we 
have that 

(3.17) 7 w / = 20u'- 10//-9 if u V^+i. and 7<Tm+1 = 20<^+i - 10/i - 19. 

5.6. The equality (|3.7[) shows that B n / u /\^ iX + C XtU >$ > A^jy)^^; item 5.5 
therefore implies 

(3.18) C XtU >$>~/u> for every u' G {o^ + 1, ... ,o>+i}. 

6.1. From Definition 13.71 it follows that A^^ = 10a — 3, so the equality (|3.7[) 
implies i?6*,4 + C^o-fl > 10c — 3. By (|3.14[) . it holds that C^ G ^ = 0, so we obtain 

(3.19) B 6 * A > 10(7 - 3. 

6.2. From Definition 13.71 it also follows that A-j^y^ — j u > for u' G {07, + 
1, . . . The equality (I3.7[) therefore implies £>7<|»,6 + Ce.u'tl > 7«'- By (|3.13l) . 
it holds that = 0, so we have 

(3.20) C 6 ,«'tt > T«' for every u' G {cr^ + 1, . . . , er^+i}. 

6.3. Again, Definition 13 . 71 implies that for every u' G {a^ + 1, . . . , <7 M+ i}, we have 
A v (u')<>,u"b = 7«' (where the function v is defined with respect to Remark I3.5|) . 
Therefore, by the equality (|3.7|) . we obtain -B„(„')o,3 + C3,«'J > 7«'- By (|3.15j) . it 
holds that B^/)^ 3 = 0, thus, 

(3.21) Cs.u'tl > lw for every u' G {cr M + 1, . . . , cr^+i}- 

6.4. We use Definition 13.71 again and obtain A v ><y^ = p v t for every v' G {<j^ — 
fj, + 1, ... , er^_|_i — (i — 1}; the number p v > is to be understood here with respect 
to (|3.1I) . namely, 

(3.22) = 2(h/ + 10/i + 1. 

Therefore, by the equality (|3.7[) . we obtain B„'(?,4+C4,8t > p^,, B„'<?,5+C5,8t > p v i. 
By (|3~T4|) . it holds that C 4i8l , = C 5 , 8 i, = 0, thus, 

(3.23) B v >Q A >Pv', B v >v£>p v > for every v' G {a^ - p,+ 1, . . . , cr^+i - /x - 1}. 

6.5. We have -^^.(a- +i)j = 10er M + 11 by Definition 13.71 Therefore, by the 
equality p.7[) , we obtain B^^ + C^r a +1)^ > lOoy, + 11. By the equality (|3.13p , it 
holds that #54,5 = 0, thus, 

(3-24) C 5X ^ +1)t > 10^ + 11. 

Step 7. 
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7.1. We have Aq^^ = 10cr M + 7 by Definition 13.71 

7.2. The conditions (|3~2H|) . (JS2H), and p~2Ij) show that C^ ((T +1)jJ > 
lOcr^ + 7 for t G {y, 3, 5, 6}; the conclusion of Step 3 that B 6 ^ £-> 0, thus we have 

( 3 - 25 ) B 6+t+ C t,(<?K + l)$ > ^6*,(a M +l)|- 

7.3. By (|37T9|) . 5 6<M > 10cr - 3. From Remark it follows that a - > 2, 
so we obtain £?6*,4 > lOo^, + 7. The equality p. 141) implies C 4 r CT +i)fl = 0, so that 

(3.26) B 6 + A + 64,(^+1)1) > Ae^^^+x)!). 

7.4. From p~TT]) it follows that, finally, 

(3-27) B 6 +j + 67,(^+1)}! = B 6 *,8 + C 8i((Tai+1)s = 00. 

7.5. The conditions (|3.25p - (|3.27[) show that the minimum in (|3.7[) for i = 6<|k, 

j = (cr^ + 1)[| is provided by r = ip (where ip stands for the index from {1,2} not 
equal to x, that is, = {1, 2} \ {x})- So, using item 7.1, we see that 

Bq+^ + C ¥ , ! ( . (j+1 )j = lOcr^ + 7; 

by #6*^ - 0, thus, 

(3-28) CV,(-,+D» = 10ct m + 7 - 

Step 5. Now we consider indexes u G {07, + 1, . . . , o>+i} and v G {<j^ — p + 
1, . . . , cr^+i — fi — 1} that satisfy the condition either v = u — fi — 1 or tJ = m — /i. 

8.1. Definition 13.71 shows that Ay^^ = 10m + 10V — 10 if u ^ cr^+i, and 

^(<T tI+1 - M -i)^, ( T M+1 tt = 20(7^+1 - 10/i - 21. 

8.2. Step 3 shows that By<y yt > and > for t G {1, . . . ,6}, the condi- 
tion p.llj) that ByV t 7 = = 00. 

8.3. The condition p.23[) shows that By^> i4 > py, 5 > Pw- Using item 8.1 
and the equality (I3.22j) . we obtain By<y 4 > Ay<y^, By<y 5 > Ay<y^. Further, the 
conditions (|3.18[) . p.20[) . and (|3.21[) imply that the elements C X) „j, C 3 ^, and Cq^ 
are all greater than or equal to 7^. Item 8.1 and the equation p,17[) show therefore 
that (7 x ,„jj, C 3 , S jj, and C 6 ^ are greater than iw^. 

8.4. Items 8.2 and 8.3 show that By<yj + Ct,™]) > Ay^y^t whenever t G 
{x, 3, 4, 5, 6, 7, 8}. The equality (|3.7[) thus shows that Byty iV + Cp,sj) = Ay^^t,, 
and, therefore, 

(3.29) By<v tV + — 10u+ 10u — 10 for u^<j^ +1 , and 

(3.30) Cy.^+ifl = 20cr AI+ i - 10/i - 21 - -B^^-p-i)^. 

8.5. By the initial setting of Step 8, v can be either u — (i or u — /1 — 1 if 
u £ {(T M + l,(T M+ i}, and v = — p + 1 if u = + 1 in ()3.29|) . Therefore, 
the conditions p.29[) together with p.28p and p.30[) form a classical linear system 
with 2(cr /i +i — er^) — 1 real unknowns and with the same number of equations; the 
system has a unique solution and it is given by By<y iV = lOv — 7, C Vt u$ = 10u — 3 
for u ^ o>+i, and 

(3.31) C^ +1 « = 10^+1 - 4- 

8.6. Finally, we note that Be^, v — by the equality p,13j) . so the equality p. 311) 
implies Bq^ j(P + C<p !(Tf4+ ijj = 10<7 M +i — 4. Now we use the equality p.7j) to derive 
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that j46<d,er„ +1 j) < lOcr^+i — 4 and obtain a contradiction with Definition 13.71 with 
respect to which Ae^ ta ^ i+1 j = lOer^+i — 3. 

The contradiction obtained shows that the initial assumption was false. There- 
fore, matrices B G S |I|x8 and C G M 8X '' 7 ' can not satisfy B ®C — A(£,7£). □ 

Now we are ready to prove the main results of the present section. 

Theorem 3.12. The k-TMF problem is NP-complete for any integer k > 8. 

Proof. By Lemmas 13.101 and 13.111 the matrix ^4 (£,7^) has factor rank at most 8 
if and only if Problem 13.21 returns a positive answer on the query (£,7£). Using 
Proposition ^. 61 we construct the matrix A(T,, 1Z, k), which has factor rank at most 
k if and only if Problem 13.21 returns a positive answer. By Lemma 13.81 the number 
of operations needed to construct A(E, 1Z, k) is polynomial in the size of the data. 
Finally, we perform an appropriate scaling and apply Lemma 12.51 thus giving a 
polynomial transformation from Problem 13.21 to fc-TMF. That fe-TMF is in NP 
follows from Theorem 12.41 □ 

Theorem 13.121 allows us to obtain a negative answer for Question 11.21 

Corollary 3.13. It is an NP-hard problem to decide whether the factor rank of a 
given tropical matrix equals 8. 

Proof. In order to check whether the factor rank is at most 8, we need to solve the 
8-TMF problem. So the result follows from Theorem l3~T2l □ 

Another result that can be obtained as a corollary of Theorem 13.121 gives an 
answer for Question 11.51 In contrast with the situation of the Gondran-Minoux 
rank, which can be computed in polynomial time if a given matrix has bounded 
tropical rank [18] . Qucstion ll.5l is answered in the negative. 

Theorem 3.14. Given a matrix A of tropical rank at most 8, it is NP-hard to 
decide whether the factor rank of A is at most 8. 

Proof. Using the definition of tropical rank explicitly, one can check in time 
0(m 9 n 9 ) whether the tropical rank of an m-by-n matrix is greater than 8. Further, 
Theorem 1.4 of [TB] shows that the factor rank of a matrix is greater than 8 if the 
tropical rank is. Thus the present theorem follows from Theorem 13. 121 □ 

4. When TMF is easy: the case k < 3 

In this section, we discuss the computational complexity of tropical matrix fac- 
torization in the case when the dimension k is small. It turns out that the TMF 
problem admits an easy solution if k is required to be small enough, indeed, the case 
of k = 1 is trivial. In |15) , the authors show that every tropical matrix whose factor 
rank is greater than 2 has a 3-by-3 submatrix with full factor rank, thus giving a 
polynomial-time algorithm for solving TMF in the case k < 2. In fact, the TMF 
problem with k < 2 can be solved by a linear-time algorithm, see [101115] . However, 
the question of whether the TMF problem can be solved in polynomial time in the 
case k = 3 remained open. The goal of this section is to give an affirmative answer 
to this question. 

The proof of the main result of this section is based on a reduction to a problem of 
linear programming. In the most general setting, the problem of linear programming 
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asks whether a given system of m linear inequalities with n real unknowns has a 
solution. The polynomial-time algorithm for this problem has been constructed by 
Khachiyan in |27j . However, the number of operations that Khachiyan's algorithm 
needs to perform depends on the magnitude of coefficients, and this number can be 
unbounded even if m and n are fixed |30j . On the other hand, it turns out that in 
our case, it is sufficient to consider the following special case of linear programming. 

Problem 4.1. Given a system of m linear inequalities involving n real unknowns 
but no more than two unknowns per inequality. 

Question: Does there exist a solution of the given system? 

In contrast with the general problem of linear programming, a strongly- 
polynomial algorithm is known for solving Problem l4.ll In other words, the number 
of operations this algorithm requires is bounded by a polynomial of m and n. 

Theorem 4.2. |30| There is an algorithm that solves Problem \4-l\ and requires at 
most O (mn 3 log rnj arithmetic operations. 

In order to prove our main result, we also need the following technical lemma. 

Lemma 4.3. Letri be the greatest element of the ith row of a matrix A G R pX9 , and 
Cj the greatest element of the jth column. If numbers Ui andvj satisfy mm{ui,Vj} — 
Aij for any pair of indexes, then either Ui = fj for every i or Vj = Cj for every j. 

Proof. If Ui' < Tii for some i', then we have Ai>j < vy for every j, thus obtaining 
a contradiction with the maximality of ry . So we see that Ui > ri for every i. 
Analogously, we can show that Vj > Cj for every j . 

Finally, if we have Ui > ri and Vj > Cj for some i and j, then A^ is greater than 
min{ri,Cj}, again contradicting the maximality of r*j and Cj. □ 

Now we will prove the main result of this section. Throughout the proof, when 
we say that a certain computation can be made in time tp(m, n), we mean that the 
computation requires at most <p(m, n) arithmetic operations in the worst case. 

Theorem 4.4. Given a matrix A G R mxn . There is an algorithm that de- 
cides whether the factor rank of A is greater than 3 and requires at most 
O ((m + n) 3 mn log(mn) + m 3 n 3 ) arithmetic operations in the worst case. 

Proof. We need to decide whether the system of equations 



has a solution B £ R mx3 , C G K 3xn . 

1. Proposition 6.1 from [15] states that A has factor rank at least 3 if and only 
if there is a full-rank 3-by-3 submatrix of A. The 3-by-3 submatrices of A can 
be treated in time 0(m 3 n 3 ), so we can assume without a loss of generality that 
A[l, 2, 3|1, 2, 3] has factor rank 3. Then we perform an appropriate scaling and 
relabeling, and can therefore assume that is the minimal element of every row 
and every column of A, and that A[l, 2, 3|1, 2, 3] equals either 



where the asterisks stand for certain (possibly different) positive numbers. 



3 



(4.1) 



min{_B ir + C T j} = A, 



(4.2) 
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2. Assume that the jth column of A consists of zero elements. In this case, we 
check whether minimi, A^, ^13} = for every i. If so, the jth column of A is a 
tropical sum of the first three columns, and we can remove the jth column without 
changing the factor rank. Otherwise, if min{^4ii, Ai2, A^} > for some i, then the 
submatrix A[l, 2, 3, i|l, 2, 3, j] equals either 



V 











/ 



V 











/ 



so the factor rank of A is at least 4. Note that zero columns can be straightforwardly 
detected in time 0(mn). We can therefore assume that no column (as well as no 
row, analogously) of A consists of zero entries. 

3. Suppose that matrices B and C indeed satisfy (|4.1[) . Then for every r G 
{1, 2, 3}, we subtract b T , the least element of the rth column of B, from every entry 
of that column and add b T to every entry of the rth row of C. The matrices B' 
and C obtained satisfy B' ®C = A, and the minimal element of every column of 
B' equals 0. From item 1 it then follows that every element of C is nonnegative as 
well. Thus we can assume without a loss of generality that 

(4.3) B„ > 0, C rj > 

for every triple of indexes i, j, r. It is thus sufficient to decide whether there exist 
matrices B G ]R mx3 and C £ R 3x ™ that satisfy flS) and 03J simultaneously. Let 
us start to describe the algorithm that makes this decision. 

4. Items 1 and 3 allow us to determine the 'zero patterns' of B[l, 2, 3|1, 2, 3] 
and C[l, 2, 3|1, 2, 3] up to relabeling the columns of B and relabeling the rows of 
C . Namely, the problem reduces to the following three cases, each of which can be 
treated by the algorithm separately: 

* * 

(4.4) £[1,2,3|1,2,3] = C[1,2,3|1,2,3]= [ * * 

■ * 



(4.5) 



5)1,2,311,2,3] = 



* * \ / * 

* * , C*[l,2,3|l,2,3] =0*0 

* * o / loo* 



(4.6) £[1,2,311,2,3] 



C[l, 2,3)1,2,3] 



If B[l, 2, 3 1 1 , 2, 3] has the form as in ()4.4)) . then, by (|4.ip . we necessarily have 
C T j = if A Tj = and C TJ > if A Tj > for every r 6 {1, 2, 3} and j € {1, . . . , n}. 
Further, if C[l, 2, 3|1, 2, 3] has the form as in ([P]) . then (H7T) implies that B iT = 
if Ai T = and Bi T > if A- lT > for every r G {1, 2, 3} and i € {1, . . . , m}. 

Now consider the case (|4.5|) , then the sign of every entry of C is determined 
uniquely by the considerations of the previous paragraph. From item 2 it follows 
that for every i, there exists j — satisfying Aj(i) > 0. If G {1,2,3}, 
then we use (|4.1|) to obtain that B iT = if C t j(q > and B iT > if C T ju) = 0. 
Otherwise, if An = A& = A^ = 0, then we have > 3, and from (|4.1|) it can be 
deduced that again B iT = if C t jm > and B iT > if C T ,j(i) = 0. 
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Thus we have determined how the 'zero patterns' of B and C look like in the 
cases (|4.4p and (|4.5p . and we have performed 0(mn) operations to do this. The 
case (|4.6p can be considered in the same way as (|4,5[) , 

5. Let Ay equal fy, the greatest element of the jth column of A. By item 4, 
we know a r(j) € {1,2,3} for which B iT ^ = 0, the equations (|4. 1 [) then imply 
that C T tj\j > f3j. Moreover, if the equations (14. 1|) and (|4.3[) are satisfied by some 
B and C, then they can not be then broken by setting C T ^j = (3j. Note that 
computations of item 5 also require at most 0{mn) operations. 

6. Thus for every j, we can assume without a loss of generality that we have 
yet computed the values C Tl u\j = 0, C T2 ^j = /3j, where Tt(J),T2(J) G {1,2,3}. 
Analogously, for every i, it be can assumed that B iiT ^ = 0, B i T ^ = jj for 
r[(i), t^i) G {1, 2, 3} (ji stands here for the greatest element of the ith row of A). 
Note that the numbers T%(j) and t[(i) are determined by the computations of item 
4, the numbers t% (j) and T 2 (i) by the computations of item 5. Note that by item 
2, f3j ^ and 7, ^ for every j and j, so that T\{j) ^ T%{j) and t{(i) ^ i~ 2 (i). We 
also denote {r 3 (j)} = {1, 2, 3} \ {ri(j), r 2 (j)} and {^(i)} = {1, 2, 3} \ {r{(i), t£(*)}. 

7. Let t±,t2,ts be different elements from {1,2,3}. By Iti,t 2 we denote the set 
of all % satisfying r{(i) = t\, t' 2 {i) — t 2 ; by J^,^ denote the set of all j satisfying 
n(j) =*i, r 2 (j) =t 2 . 

Consider the subsystem of the system (|4.1[) formed by the equations with i G I\ : 2, 
j £ ^3,2 • We can now apply Lemma 14.31 and reduce our problem to the two cases, 
each of which can be treated separately. In one of these cases, the value of B^ is 
determined for every i G ^1,2] in the other of them, C\j is determined for every 

3 G Ja.2- 

We repeat our considerations with arbitrary Z tljt2 and Jt 3 ,t 2 instead of I\,2 and 
Jz,2, the algorithm needs then consider 2 6 = 64 cases separately. In what follows, 
we can therefore assume that, if 73 (i) = T 2 {j) and r{(z) ^ then either B i T >^ 

or C T3 (j\j is yet known. Finally, we note that the computations of item 7 can also 
be performed in time 0(mn). 

8. Now we will prove that for every i and j, the equation (|4.1|) either gives 
a contradiction or is equivalent to a linear inequality or equation with at most 2 
unknowns. This will reduce our problem to Problem l4.H and then we will complete 
our proof by applying Theorem 14.21 

Indeed, let i G {l,...,m}, j G {l,...,n}. If the entry B ijT i^ is known, 
then (|4.ip is either a contradiction or can be rewritten as either i? ijT3 y) +C T3 (j), j — 
A^ or _B iT3 (j) + Ct 3 (j)j > Ay. The case when the entry C T3 y)j is known can be 
treated in the same way. There are the four possibilities in the case when both 
Bi^d) and C , T3 ( J )j- are unknown. 

Case A. Suppose r^(i) = r 3 (j). In this case, if (|4.1[) is not a contradiction, then 
it can be rewritten as either B i T ^ + C T3 (j),j = Ay or B i<T >^ + C T3 ^j > Aij. 

Case B. Let T^(i) = t^O)- By item 4, we have B i T ^ > 0, thus B iT ^ + 
C T ^(i)j- > /3j > Ay. If (|4.ip is not a contradiction, then it can be therefore rewritten 
as either B hT3(j) + C T3(j)tJ = Ay or B iMj) + C Ts(j) j > Ay. 

Case C. Let r^i) = ^(j). By item 4, we have B iT3 ^ > 0, thus Bj jT3 y) + 
> 7* — ■ If <|4.1|) is not a contradiction, then it can be therefore rewritten 
as either Bi, T |(i) + = or s i,T^(i) + > A/- 

Case D. The possibility ^(i) = T2(j), t{(i) ^ Ti(j) is not indeed realizable 
thanks to item 7. 
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The cases A-D cover all the possibilities, so that for every indexes i and j, the 
equation (|4.1j) cither gives a contradiction or is equivalent to a linear inequality or 
equation with at most 2 unknowns. What we need is to decide whether the sys- 
tems (|4.ip and (|4.3|) have a simultaneous solution, and this can be now done in time 
O ((m + n) 3 mnlog(mn)) by applying the algorithm mentioned in Theorem l4.2l □ 

Theorem 14. 41 shows that the TMF problem admits a polynomial-time solution in 
the case k = 3. We can now obtain a positive answer for Question II .31 

Theorem 4.5. There is a polynomial-time algorithm that decides whether a given 
tropical matrix has factor rank 3. 

Proof. Checking whether the factor rank is at least 3 requires at most 0(m 3 n 3 ) 
operations because it is sufficient to check whether some 3-by-3 submatrix has 
factor rank 3, see Proposition 6.1 from [15]. Checking whether the factor rank is 
at most 3 can be done in polynomial time as Theorem 14.41 shows. □ 

We also prove the result on the complexity of the 3-TMF problem. 

Theorem 4.6. The 3- TMF problem belongs to the class P. 

Proof. Solving 3-TMF means deciding whether a given matrix has factor rank at 
most 3. So the result follows from Theorem 14.41 □ 

The fc-TMF problem, which belongs to NP by Theorem l2.4[ is still not known to 
be either in P or NP-complete if k £ {4, 5, 6, 7}. For these values of k, the question 
of whether /c-TMF admits a polynomial-time solution seems to be an interesting 
topic for future research. 

5. Factor rank and minors of fixed size 

In this section, we discuss a problem that is related to determining computational 
complexity of tropical matrix factorization. The problem we are interested in, 
Question 11.41 of the present paper, has been formulated by Develin, Santos, and 
Sturmfels in [15] . The question asks to find some N = N(k) for which every matrix 
with rank greater than k necessarily contains a submatrix whose rank is greater 
than k and whose size is at most N x N. For instance, the basic result of classical 
linear algebra states that we can set Nr(k) = k + 1 for usual rank function of 
matrices over a field. Also, it follows directly from the definition that for tropical 
rank of a tropical matrix, one can set Nt r (k) = k + 1. The same question has also 
been studied for the Gondran-Minoux rank functions of tropical matrices, and it 
has been proved that setting NcM(k) = 2 fe+1 is sufficient in that case [15] . 

For the factor rank, it is trivial that we can set N(0) = 1, and it is easy to 
check that setting N(l) = 2 suffices. The fact that N(2) can be set to equal 3 has 
been proved in [15] . It was not known, however, whether there exists an N(k) if 
k > 3. If there were such N(k), it would be bounded below by a function having 
an exponential growth as k approaches infinity, see Proposition 2.2 of |15j . In this 
section, we construct an example showing that indeed no such N(k) exists for k > 4. 

We fix an arbitrary integer v > 1 and will use its designation throughout this 
section. Our goal is to construct a matrix C with the following properties: (i) every 
minor of C of size at most v x v has factor rank at most 4, (ii) the factor rank of 
C itself is greater than 4. Now we present an example of such a matrix. 
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Example 5.1. Denote an arbitrary number from {1,...,^ + 1} by u, from 
{1, . . .,v + 2} by v. The entries of the tropical matrix C = C{v) £ r("+6)x(»+6) 
are defined as follows: 

(1) C uv = min {10u + lOv — 10, 20u + 1, 20v — 9} if{u, v) equals neither {v, v+1) 
nor (v + 1, v + 1); 

(2) C VjU+ i = 20is - 1; C v+ljU+ i = 20is + 11; 

(3) C M , y +3 = 0, C u ^+4 — C MjI y+5 = C Mj y+6 — 2; 

(4) Cv+2,u = 10m — 3, C u +3^ u = 11, CV+4,u = C y +5 !tl = 2, = 0; 

(5) the matrix C[v + 2, . . . , v + 6|f + 2, . . . , v + 6] equals 

( 2 \ 

2 2 

2 2 2 2 . 

2 2 2 

\ 2 2 2 / 

Let us prove that the factor rank of the matrix introduced is at least 5. 

Theorem 5.2. The factor rank of the matrix C(y) is greater than 4. 

Proof. Assume the converse. Definition 11.11 then implies that C = A ® B for some 
A e R(^+ 6 ) x4 and B € that is, 

4 

(5.1) min{A sg + B qt } = C st 

q=l 

for every s, t 6 {1, . . . , v + 6}. 

1. Note that the product A® B will remain unchanged if we add some number 
r£tto every entry of the qih column of A and — r to every entry of the gth row of 
B. So we can assume without a loss of generality that for every q € {1, 2, 3, 4}, the 
qth row of B has as its minimal element. In this case, the matrix A also consists 
of nonnegative elements, because from Example 15.11 it follows that the elements of 
C = A <E> B are nonnegative. 

2. Let us denote by A' the 4-by-4 matrix formed by the rows of A with labels 
v + 1, v + 4, v + 5, v + 6, by B' the 4-by-4 matrix formed by the columns of B with 
numbers v + 3, v + 4, v + 5, v + 6. From the definition of the matrix C, which is 
equal to A (g) B, it follows that the product A' (g) B' has the form 



(5.2) 



( 


* 


* 









* 















* 


* 


o 1 



where the asterisks stand for certain (possibly different) positive elements. Item 1 
shows that A' and B' are nonnegative, so that the matrices A' and B' have the 
form (|5.2p up to permutations of rows and columns. Note that performing the 
same permutation on column indexes of A and row indexes of B will not change 
the product A®B. So we can assume without a loss of generality that the matrices 
A' and B' have exactly the form (|5.2p . 

3. From item 1 it follows that the elements of A and B are nonnegative, from item 
2 that A u+5tq — if and only if q — 3. So from (|5.ip it follows that B 3t = if and 
only if Cy+5.t = 0, where t is an arbitrary index from {1, . . . , v + 6}. Analogously, 
item 2 implies that A V +§ A — if and only if q = 4, so that the equation (|5.ip shows 
that Bn = if and only if C u+ e tt — 0. 



22 



YAROSLAV SHITOV 



4. Item 2 also implies that = if and only if q = 1; Bq jV +4, = if and 
only if q — 2; B qM +$ — if and only if g = 3. Thus the equation (|5.1[) shows that 
A s \ = if and only if = 0; A S 2 — if and only if C Si „+4 = 0; A s a = if and 
only if = 0; here s runs over {1, . . . , v + 6}. 

5. In particular, items 3 and 4 imply that the entries A v +i,i, A„+3 t x, A v+ 2.2i 
A u+? ,^, i?3^+2, B i>1/+ x, B i}U+2 are all equal to 0. 

6. From (|5.1|) it also follows that min{A s3 + B 3yl/+ 2, A s4 + B itU+ 2} > C StU+ 2 
for every s G {1, + 6}. Item 5 implies that -63,^+2 = B^ u+ 2 = 0, so we 
have min{A S 3, A s ±\ > C s _ u +2- Applying the definition of the matrix C, we obtain 
Ai 3 > 20i + 1, A iA > 20i + 1 for every i G {1, . . . , v). 

7. Analogously, from item 5 it follows that A„+i t i ~ 0, so (|5. 1[) implies that 
Bu > Cw+i,t for every t G {1, . . . , v + 6}. For j 6 {1, . . . , v + 1}, we therefore 
obtain By > 20 j - 9. 

8. Again, item 5 implies B^+i = 0, so (15. ip shows that A u+ 2a > CV+2,1,,+1 = 
10^ + 7. Further, by item 5, A u+3 ^ — A u+3y3 = 0, so (|5.ip implies min{_Bn, -B31} > 
C u+3 ,i, that is, Bu > 11, B 31 > 11. 

9. Item 8 shows that the number A v+ 2, q + B q \ is greater than 7 unless q = 2. 
The definition of the matrix C shows that = 7, so from (|5.1j) it follows that 
A u+ 2,2 + B21 = 7. Item 5 implies that A^+2.2 = 0, and we therefore obtain B21 = 7. 

10. Now let i £ {1, . . . , v} and either j = i or j = i + 1. The definition of the 
matrix C implies CVj < lOi + 10 j — 10. From items 6 and 7 it then follows that 
the numbers A^, An, and -By are always greater than CV,. The equation (|5.ip 
therefore implies A^ + B2j = Cij. We thus obtain the following system of 2v + 1 
linear (in the classical sense) equations with 2^ + 1 real unknowns: B21 — 7 and 

(5.3) Ai2 + B 2 i = Cu, A i2 + B 2 , l+ \ = Ci t i+x, for every i G {1, . . . , v}. 

Solving the system, we obtain that, in particular, B2.U+1 = 10f + 6. From the 
equation (|5.ip it thus follows that C v+ 2m+i < ^+2,2 + B 2 ^+i- Item 5 shows that 
Av+2,2 = 0, so that C u +2,v+i < 10f + 6. However, the definition of C shows that 
C v +2.v+\ — 10f + 7. The contradiction obtained shows that our initial assumption 
is not true, so the result follows. □ 

Let us show that every v-\sy-v submatrix of C(y) has factor rank at most 4. 

Theorem 5.3. Let fi 6 {1, . . . , v + 1}- The matrix C that is obtained from C{v) 
by removing the /ith column has factor rank at most 4. 

Proof. We need to provide matrices A G R( ,y + 6 ) x4 and B G R 4x ^+ 5 ) satisfying 
A(E)B = C . Denoting an arbitrary number from {1, . . . , z/+l} by i, from {1, . . . , v+ 
1} \ {fi} by j, we can construct these matrices as follows. 

We set An = and A& = An = 20i + 1 for every i, Ai2 — lOi — 7 if i < fi, 
A i2 = Wi - 8 if (j, < i < v + 1, and A v+1 .2 = 2Qv + 21. Set also B xj = 20j - 9 and 
B i3 = for every j, B 3j = Wj - 3 if j ^ 1, B aj = 11 if j = 1. Define B 2] = lOj - 3 
if either j < fi or j = v + 1, set B2j = 10 j — 2 if /j, < j < v + 1. 

Finally, we define the matrices 



/ 

11 

2 2 

2 2 



10^ + 7 



11 

2 
2 




\ 



and 



/ 20^ + 20 2 2 

20^ + 20 2 2 

2 2 

\ 2 2 2 



/ 



2 \ 



2 
2 



\ 2 2 2 
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to be, respectively, the matrix formed by the rows of A with numbers u+2, . . . , v+6 
and the matrix formed by the columns of B with numbers v + 2, . . . ,v + 6. The 
product A £g> B can now be straightforwardly checked to equal the matrix C. □ 

Note that the factor rank of a submatrix is always less than or equal to the factor 
rank of whole matrix, Theorem 15.31 thus implies that every submatrix of C (v) of 
size at most v-by-v has factor rank at most 4. Therefore, the matrix C(v) gives 
a negative answer for Question 11.41 if k — 4. We are now ready to prove the main 
result of this section, and answer Question ll.4l for every k > 4. 

Theorem 5.4. Let k and v be integers and k > 4. There exists a tropical matrix 
C(k,v) with the following properties: (i) the factor rank of C(k,v) is greater than 
k; (ii) every minor ofC(k,v) of size at most v-by-v has factor rank at most k. 

Proof. Apply Proposition ^. 6l and then Lemma 12.51 to the matrix from Example l5.ll 

□ 
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